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OPERATIONAL PROCEDURE FOR COMPUTER PROGRAM FOR DESIGN-POINT 
CHARACTERISTICS OF A GAS GENERATOR OR A TURBOJET LIFT ENGINE 

FOR V/STOL APPLICATIONS 
by Richard P. Krebs 
Lewis Research Center 

SUMMARY 

The computer program described in this report calculates the design-point charac- 
teristics of a gas turbine engine which can be used as a gas generator or a turbojet lift 
engine. Intended applications are in the field of propulsion for V/STOL aircraft. 

The program computes the dimensions and mass, as well as the thermodynamic 
performance of a model gas turbine engine. Physical and thermodynamic character- 
istics of the engine components are also given. The report describes the engine model 
used in making the calculations. 

Details on the required input parameters are given, together with a description of 
the preparation of the input data cards. The option for selecting either a gas generator 
or a lift engine design is discussed, and methods for running several designs with a 
single submission of the program deck are illustrated. 

The program was written in FORTRAN IV language. Provision has been made so 
that the program accepts input values in either SI Units or U. S. Customary Units. Each 
design-point calculation requires less than 0. 5 second of 7094 computer time for execu- 
tion. 


INTRODUCTION 


The Lewis Research Center is interested in propulsion systems for direct lift in 
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attractive candidates for these applications. The lift fan category has been broken down 
into the integral lift fan system in which the drive engine is directly connected to the 
fan, and the remote -drive lift fan system in which the lift fan is driven by a turbine 
mounted at the tip of the fan blades (e. g. , ref. 1). The latter has been further divided 
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into systems where the tip turbine is driven by the exhaust from a gas turbine engine 
(gas generator), and systems where the tip turbine is driven by compressed air from 
an air generator. 

As a part of an analytical research effort in the field of direct-lift propulsion sys- 
tems, a number of electronic digital computer programs have been written to determine 
the design-point thermodynamic characteristics, geometry, and mass of these systems. 
Two of these computer programs have already been described and reported in refer- 
ences 2 and 3. The computer program of reference 2 provides a preliminary design 
and analysis tool for an entire tip-turbine -driven lift fan assembly. This program is 
useful for the initial sizing of the fan and turbine. It is also adaptable to parametric 
studies of the effect of changes in the principal design variables of both the fan and 
turbine on the design-point characteristics of the fan assembly. 

Reference 3 describes a computer program which calculates the design-point char- 
acteristics of a compressed air generator such as would be used to power a tip-turbine- 
driven lift fan. The program computes the dimensions and mass, as well as the thermo- 
dynamic performance, of a model air generator configuration consisting of a low spool 
with its air compressor and drive turbine and a high spool which serves as a gas gener- 
ator for the drive turbine on the low spool. Physical and thermodynamic characteristics 
of the air generator components are also given. 

The present report deals with two other components for V/STOL propulsion systems. 
A single computer program, described herein, contains an option whereby the design- 
point characteristics of either a single spool gas generator or a gas turbine lift engine 
can be calculated. When the gas generator option is used, the size is established by 
the prescribed gas flow. In the other option, the size of the lift engine is determined 
from a prescribed jet thrust. The analysis, approach, and procedure are similar to 
those used in the air generator program described in reference 3. 

This report was compiled to furnish the descriptions and instructions necessary 
for running the computer program for either option indicated in the preceding paragraph. 
It was assumed that the user of this report is familiar with digital computer program- 
ming and is knowledgeable concerning the parameters used in describing gas turbine 
engines. A complete description of the input parameters is included, together with in- 
structions as to how the input data are preared. A typical computer output page is in- 
cluded, and the meaning of each output parameter is given. The FORTRAN program 
statements are included for those who wish to know more about the program, or for 
those who may wish to change it. 
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APPROACH 
Engine Model 


The component arrangement and station locations for the gas turbine configuration 
upon which the development of the computer program is based are shown in figure 1. 
Aside, from the inlet and exhaust sections, the gas generator and jet lift engines are 
assumed to have the same component arrangements in the analysis used for the com- 
puter program. From the inlet, air enters the compressor where it is compressed by 
a multistage axial -flow unit. The largest part of the high-pressure air upon leaving the 
compressor enters the combustion chamber where it mixes with the fuel, and the mix- 
ture is burned. If the temperature of the gas leaving the combustor is sufficiently high 
so that the turbine needs to be air cooled, the cooling air may be taken from the com- 
pressor discharge. Air may also be bled from the compressor for aircraft control and 
other purposes. The hot gas from the combustor passes through the turbine which may 
have one or more stages. The bare (core) engine, then, from the compressor inlet face 
to the turbine discharge is the same in concept for both the gas generator and the jet 
lift engine. 

Because the engine inlet is so intimately related to the engine installation, and 
because each inlet will be different for each installation, no attempt was made to define 
the inlet for either type of engine. In general, it is probable that the inlet for a gas 
generator would be long enough to accomodate some form of acoustic treatment in order 
to reduce the noise generated by the compressor. This assumption is probably true 
whether the gas generator is forward or rearward facing. On the other hand, for a lift 
engine, minimal length is of extreme importance, and the compressor noise would be 
masked by the jet noise so that inlet acoustic treatment may not be necessary. Accord- 
ingly, the inlet of the lift engine would probably be as short as possible, commensurate 
with introducing the flow to the compressor with an acceptable maximum flow distortion 
and energy loss under the most adverse flight conditions. 

On the exhaust end, it is somewhat easier to define a minimal hardware require- 
ment for both the gas generator and the lift engine. In the case of the gas generator, a 
short duct is required to diffuse the turbine flow velocity down to a velocity such as could 
be used in a gas transfer line between the gas generator and the tip-mounted turbine 
driving the lift fan. This diffusion should take place with a minimal loss to conserve 
the pressure energy to drive the fan. 

The model of the exhaust duct assumed for the gas generator is shown in figure 2. 
For simplicity in the weight and length calculation, a fixed duct geometry was assumed 
that was believed to be representative of actual configurations. The following quantities 
were defined in terms of the rotor tip diameter at the turbine discharge D^.: 
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Axial length: = 0. 75 D^. 

Exit diameter: = 0. 75 D^. 

Turbine exit hub-tip ratio: D^/D^. = 0. 75 
The weight of the outer shell plus the inner cone was then approximated by 

w d - K d D , 2 


0 ) 


For representative input values such as a material density of 8300 kilograms per cubic 
meter (520 lbm/cu ft) and an effective wall thickness of 0. 85 millimeter (0. 033 in. ), it 
was calculated that = 23 (K^ = 4. 7 for D^. in ft, in lbm). 

On the lift engine, the exhaust section takes on the function of a nozzle and speeds 
up the flow so that the static pressure at the end of the nozzle is ambient or near ambient 
pressure. The acceleration can be accomplished in a much shorter length than can the 
diffusion. Again, for simplicity, a representative nozzle configuration was taken such 
that 

L n = 0.25 D t 

D = 0. 80 D. 
n t 

No central cone was considered for the nozzle. The weight of the nozzle was estimated 
as 


w „ = K D, 2 
n n t 

where K n = 5.5 (K^ = 1.1 for D^. in ft, W R in lbm). 


( 2 ) 


Computer Program 

The computer program - ONE SPOOL GAS GENERATOR/LIFT ENGINE - LIFT 
COMPONENTS - provides a design point configuration for the gas generator model shown 
in figure 1 or for a lift engine with a similar component arrangement. The gas generator 
is sized for a prescribed gas flow, while the lift engine is sized to produce a prescribed 
thrust. The overall thermodynamic performance is described for either engine along 
with the overall dimensions and total mass. Thermodynamic performance, size, and 
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mass are also calculated for the principal components. The length and mass calcula- 
tions, exclusive of the inlet and exhaust sections, are based on the lift engine compo- 
nent mass correlations presented in reference 4. 

There is a second form of the computer program described in this report - a so- 
called "cruise" version, which reflects the design changes representative of a contin- 
uously operating engine. Continuous operation would be necessary if the gas generator 
were to be used directly for cruise thrust or to supply the gas for the operation of a 
cruise fan. Component and overall masses in the "cruise" form of the program are 
calculated for the bare engine from the cruise engine component mass correlations pre- 
sented in reference 4 and yield masses that are greater than those for the lift version. 
Thus, by using the two forms of the program, it is possible to cover the range of gas 
generator component and overall masses likely to be encountered in realistic designs 
of gas generators for commercial V/STOL transport applications. In fact, a compari- 
son of calculated masses as contained in the program for two recently developed gas 
generator engines showed that the real engine mass corresponded to values of calculated 
mass at around 50 to 60 percent of the difference between the lift and cruise mass deter- 
minations. Because there would be no interest in a continuously operating and relatively 
heavy lift engine, no provision was made in the cruise version of the program to calcu- 
late the performance of a turbojet lift engine. 

The remainder of the report emphasizes the ONE SPOOL GAS GENERATOR/LIFT 
ENGINE - LIFT COMPONENTS form of the program. Any significant differences which 
appear in the cruise form are specifically indicated. 

The program was written in FORTRAN IV language for use on an IBM 7094, Model 2, 
computer. With modifications this program can be used on all machines that have a 
FORTRAN compiler. The program was developed in U. S. Customary Units, but it will 
perform the calculations for either SI inputs or U. S. customary inputs. Each pass 
through the program requires less than 0. 5 second on a 7094 computer. 

FORTRAN listings for the first form of the computer program (lift engine option, 
lift engine mass correlation) and the subroutines are shown in figure 3. 


INPUT PARAMETERS 

This computer program for the design-point characteristics of a gas generator 
engine or a turbojet lift engine has considerable inherent flexibility. Some idea of the 
flexibility can be achieved from the fact that no less than 26 independent input param- 
eters may be specified for any one engine design. In addition, values of four program 
control parameters must be supplied. The function of these four control parameters is 
discussed in the DATA INPUT CARDS section. In the following paragraphs the signifi- 
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cance of each input parameter is discussed. The symbols used in the FORTRAN lan- 
guage of the computer program are also indicated by capital letters. An attempt has 
been made to use symbols which are descriptive of the property, component, and engine 
station as indicated in figure 1. The program accepts either SI Units or U. S. Customary 
Units. Dimensions for the parameters in both systems of units are given in the DATA 
INPUT CARDS section. 


Ambient Conditions and Engine Inlet 

The ambient conditions of pressure PO and temperature TO must be specified. 
Because this is a design-point program for a powerplant that is to be used primarily for 
lift, it is assumed in this program that the total pressure and temperature at the inlet 
correspond to ambient static conditons (i. e. , the engine is designed for the takeoff con- 
dition). 

The performance of the engine inlet is described by a single parameter, the total 
pressure ratio PI2P1 across the inlet. This pressure ratio is very close to one for a 
short, clean inlet, but it may be considerably less than that if the inlet flow path is 
tortuous or filled with acoustic absorption devices. The inlet flow is assumed to be 
adiabatic. 


Compressor 

Air passes from the inlet into the compressor which has an overall pressure ratio 
equal to PC2PC1. This compression takes place in SNC stages with the corrected tip 
speed of the initial stage equal to UTIPCC. The compressor efficiency ETAC, which 
may be given in either the adiabatic or polytropic form, should reflect the dependence 
on stage pressure ratio. Relations between corrected tip speed, overall pressure ratio, 
number of stages, aerodynamic loading, and stage pressure ratio are illustrated in 
reference 4. The flow path for the compressor, whether constant hub, constant mean, 
or constant tip, is determined by the value assigned to JCGEOM. 

The tip diameter of the compressor is set by the average axial inlet Mach number 
AMC1, the inlet hub-tip ratio DHDTC1 for the first rotor row, and the airflow rate. 
Large airflow rates per unit flow area are advantageous for small engine size and 
weight. Large specific airflows are the result of high axial inlet Mach numbers and low 
hub-tip ratios. However, the aerodynamics of the velocity diagram generally limits 
these values to something less than 0. 6 and greater than 0. 4, respectively. The air- 
flow rate is determined from the gas flow rate in the case of the gas generator or the 
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required jet thrust in the case of the turbojet lift engine. Overall diffusion through the 
compressor can be regulated by the selection of the axial velocity ratio across the com- 
pressor VC2VC1. The reduction of the axial velocity through the compressor compen- 
sates for the increase in density of the air so that reasonable blade heights exist at the 
compressor discharge. The diffusion also provides for a decrease of velocity into the 
combustor. 

Two difference quantities of air bleed from the discharge of the compressor may be 
specified. One of them is used to cool the turbine and is discussed in the Cooling Air- 
flow section. The other, a so-called "user” bleed (BUSER), is available for aircraft 
control and other purposes. This user bleed is usually limited to a few percent of the 
compressor flow, but it may be up to around 15 percent in some applications. 


Combustor 

Only three input parameters are required to establish the performance and geom- 
etry of the combustor. Although the thermodynamic properties of the products of com- 
bustion are based on a fuel with a hydrogen-carbon ratio of 2, the program accomodates 
fuels with different heating values. The heating value HF is one of the independent 
combustor parameters. For JP fuel, the heating value is 42 800 kilojoules per kilo- 
gram (18 400 Btu/lb). Another independent parameter is the combustor efficiency 
ETAB. 

A third input is the combustor pressure loss PB2PB1 expressed as the overall total 
pressure ratio across the combustor. This value should generally be a function of the 
combustor length to height ratio (ref. 4). The ratio of combustor length to height and 
combustor reference velocity are fixed within the program. In the lift version of the 
program, these values are 2 and approximately 24 meters per second (80 ft/sec), 
respectively. In the cruise version, the values are 3 and approximately 18 meters per 
second (60 ft/sec), respectively, reflecting a relaxation in the severity of the combustor 
design in this latter type of application. 


Turbine 

The number of stages SNT required for the turbine depends on the compressor pres 
sure ratio and the desired turbine efficiency. However, one or two stages are satis- 
factory for most gas generator or jet lift engines. The parameter ALPHAT represents 
the angle of the flow coming outof the turbine stator as measured from the axis of 
rotation. The turbine loss coefficient AKCT sets the level of loss and, therefore, 
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efficiency in the turbine. A nominal value for the model loss relations in this program 
is between 0. 35 and 0. 40 which results in adiabatic efficiencies in the neighborhood of 
0. 88 to 0. 92 over the range of the speed-work parameter encountered in gas generators 
or lift engines. Increasing the value of AKCT decreases the turbine efficiency. 

The turbine performance is predicated for a symmetrical velocity diagram in all 
stages. Turbine efficiency is determined from this diagram and an analysis of stage 
efficiency similar to the approach used in references 5 and 6. Each stage uses the same 
constant mean diameter. The value of the mean diameter relative to the tip diameter 
at the inlet to the compressor can be set by the parameter DTMDC1, which is the ratio 
of the turbine mean diameter to the compressor inlet tip diameter. The value of 
DTMDC1 is usually less than one. A small value of DTMDC1 reduces the engine weight, 
but it also reduces the efficiency (through the speed-work parameter) and the hub-tip 
diameter ratio at the turbine discharge. The minimum value of DTMDC1 is that value 
which yields an acceptable value of turbine outlet hub -tip ratio. 

One of the most significant parameters in the performance of the gas generator or 
jet lift engine is the temperature at the inlet to the turbine stator TT1. Values used 
for this temperature should reflect allowable blade stress limits and should influence 
the amount of cooling airflow prescribed. 


Cooling Airflow 

The cooling airflow for the turbine PCA is expressed as a fraction of the compres- 
sor airflow less the amount of user bleed. The cooling air can be expressed independ- 
ently by setting the parameter PCA equal to the desired fractional value and by setting 
the engine application parameter KIND = 0. 

Two schedules of cooling air with turbine inlet temperature are also built into the 
program. The first schedule is intended to represent a cooling air requirement for an 
engine used in a lift application where the time of operation during a cycle would be 
relatively short: 


PCA = 0. 00011 *TT1 - 0. 242 (3) 

The second schedule yields a cooling air requirement typical of continuous operation at 
the assigned temperature TT1: 

PCA = 0. 00015*TT1 - 0. 297 (4) 

These two schedules are activated by setting KIND = 50 or KIND = 100, respectively. 
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Discharge Duct or Nozzle 


One of the chief differences between the two options, gas generator or turbojet lift 
engine, is the treatment of the exhaust gases from the turbine. For the gas generator 
a short duct is assumed between the turbine and the pipe directing the gas to a tip turbine 
(fig. 2). The velocity of the flow at the discharge of this short duct is established 
through the duct exit Mach number AMD2. A value of AMD2 = 0.3 represents a reason- 
able compromise between high frictional pressure drop and small duct diameter. The 
flow rate out of the gas generator, the principal parameter required in designing the gas 
generator, is WG. A total pressure ratio across the duct is included to account for any 
pressure losses associated with a particular discharge duct geometry or design. This 
ratio is given by PD2PD1. 

When the option to design a jet lift engine is used, the emphasis is on the minimum 
permissible amount of hardware downstream of the turbine. The duct is shortened, and 
it takes on the dual role of duct and nozzle. Assumptions on the size and weight of this 
duct nozzle were indicated earlier in the discussion of the engine model. 

In place of the gas flow used to size the gas generator, the jet thrust FJS is the 
sizing parameter for the jet lift engine. The duct pressure ratio PD2PD1 is maintained, 
but it will generally have a value close to 1. 0. A nozzle velocity or thrust coefficient 
is also included and is designated as CFJ. 

DATA INPUT CARDS 

The 26 independent parameters and the four program control parameters required 
as input for both the lift and cruise form of the computer program are entered on five 
data input cards. An additional card is also required which serves as an identification 
card. Although these six cards are required for a single engine analysis, other engines 
can be analyzed in a single submission of the program deck by adding one or more data 
cards in the manner described in the section Multiple Cases. The program control 
parameters N and NN direct the multiple case operation as will be explained later. The 
control parameter UNITS determines the type of units to be used in the input and output. 
The fourth control parameter OPTION determines whether the lift version of the pro- 
gram is used to design a gas generator or a jet lift engine. This parameter and function 
are not included in the cruise version of the program. 


Single Case 

In the discussion which follows, three columns are used to represent the informa- 
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tion. The first column is the name of the input variable, or parameter, as it appears 
on the computer printout sheet (see fig. 4). The second column contains the FORTRAN 
language symbol for the variable which is used in the computer program and which is 
referenced in the section INPUT PARAMETERS. The last column contains a descrip- 
tion of the variable and the units used in this program. 

The first data card sets the flow path for the compressor and selects the cooling 
air schedule to be used. The second card identifies the gas generator or lift engine. 

All the remaining data are entered on the last four cards, which use a 8F10.0 format. 

In general, the arrangement of the parameters on these last four cards is in the ascend- 
ing order of frequency of change. This arrangement can be used to advantage when 
several analyses are being run with one program submission (see the section Multiple 
Cases). 

First card . - There are four fixed-point variables used as input parameters in this 
program, and they appear on the first data card. The format for the first card is 415. 


COMP FLOW PATH 


ENGINE APPLICATION 


JCGEOM Sets the geometry of the compressor. If 

JCGEOM = 1, the compressor has a constant 
hub diameter. If JCGEOM = 2, the com- 
pressor has a constant mean diameter. If 
JCGEOM = 3, the compressor is of constant 
tip design. 

KIND This parameter selects the cooling air 

schedule used on the turbine. If KIND = 100, 
the cooling air is scheduled with a turbine 
inlet temperature to represent continuous 
operation at the assigned turbine inlet tem- 
perature. If KIND = 50, the schedule of 
cooling aii is that which might be used in a 
lift application where the time of operation 
between shutdowns would be relatively short. 
If KIND = 0, the program user may specify 
any amount of cooling air PCA as a fraction 
of the compressor flow less the user bleed. 


INPUT CARD INDEX N Set equal to 1 for single case. 

INPUT CARD INDEX NN Set equal to 1 for single case. 
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Second card . - The content of this card, which uses a 12A6 format, can be used to 
identify the particular engine case calculated. This identification is printed out on the 
third line of the computer printout. For example, on the printout shown in figure 4(a), 
the gas generator was identified as TEST CASE - CUSTOMARY UNITS. 

Third card. - Input variables which are seldom changed appear on the third card. 


TURB NOZZLE ANGLE 

ALPHAT 

Outlet flow angle, measured from the axial 
direction, for the turbine stator, in radians 
(deg). 

DUCT EXIT MACH NO 

AMD2 

Mach number which sizes the outlet of the 
gas delivery duct for a gas generator. 

PCA 

PCA 

Turbine cooling air expressed as a fraction 
of the compressor airflow less the user 
bleed. PCA is ignored unless KIND = 0 
(see first card). 

INLET HUB -TIP RATIO 

DHDTC1 

Ratio of hub diameter to tip diameter at the 
inlet of the compressor. 

TURB LOSS COEF 

AKCT 

Coefficient which sets the level of losses in 
the turbine. Nominal value, 0. 35 to 0. 40. 

FUEL HEATING VALUE 

HF 

Heating value of the fuel, in kilojoules per 
kilogram (Btu/lb). 

INLET RECOVERY 

PI2P1 

Total-pressure ratio across the inlet ahead 
of the compressor. 

INLET AXIAL MACH NO 

AMC1 

Compressor inlet Mach number. 

Fourth card. - The fourth card contains only five infrequently changed parameters. 

AX AIL VELOCITY RATIO 

VC2VC1 

Outlet axial velocity divided by inlet axial 
velocity for the compressor. 

COMBUSTOR EFFICIENCY 

ETAB 

Efficiency of combustion. 

TURB STRAIGHT VANES 

AKUT 

Index for vane or no -vane downstream of the 
turbine. If AKUT = 0, there is no vane. If 
an impulse straightening vane is desired, 
set 0 < AKUT <1.9. The loss across the 
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UNITS UNITS 

USER BLEED BUSER 

Fifth card . - 
on the fifth card. 

COMP EFFICIENCY ETAC 


FLARE RATIO (DTM/DC1) DTMDC1 
CORR TIP SPEED UTIPCC 


AMBIENT PRESSURE PO 

AMBIENT TEMP TO 

COMBUSTOR PRES RATIO PB2PB1 
NO. COMP. STAGES SNC 


vane is AKUT times one-half the loss in a 
turbine stator. 

If UNITS = 0, all quantities are in U. S. 
Customary Units. If UNITS ^ 0, SI Units 
are used. 

Bleed, expressed as a fraction of compres- 
sor flow, from peak cycle pressure for use 
outside gas generator. 


Compressor efficiency. If a positive value 
is used, the program treats it as an adia- 
batic efficiency. If the value is preceded by 
a minus sign, the efficiency is considered 
polytropic. 

The turbine is positioned radially by the 
selection of the ratio of the turbine mean 
diameter to the compressor diameter at the 
inlet tip. 

The tip speed at the inlet of the compressor, 
in meters per second (ft/sec), divided by 
the square root of the ratio compressor 
inlet total temperature to standard temper- 
ature. 

Pressure, in kilonewtons per square meter 
(lb/sq ft), for design altitude. 

Temperature, in K (°R), for design tem- 
perature situation, independent of pressure 
altitude. 

Total-pressure ratio across combustor. 

Number of stages required to produce the 
overall pressure ratio on the compressor,, 


Input variables which may be changed more or less frequently appear 
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DUCT PRES RATIO 


PD2PD1 


Total -pressure ratio across discharge duct. 


Sixth card . - What were estimated to be the most important, and, therefore, most 
frequently changed variables are read from the sixth card. 


GAS FLOW 

WG 

Flow from gas generator discharge duct, in 
kilograms per second (lb/sec). Can be 
omitted when lift engine option is used. 

COMP PRES RATIO 

PC2PC1 

Overall pressure ratio across the compres- 
sor. 

TURBINE INLET TEMP 

TT1 

Maximum cycle temperature in the engine, 
in K (°R). 

TURBINE STAGES 

SNT 

Number of stages on the turbine. 

OPTION 

OPTION 

Value of parameter selects type of engine 
to be designed. If OPTION = 1, gas gen- 
erator performance will be calculated; if 
OPTION = 2, turbojet lift engine perform- 
ance will be calculated. This option not 
available in cruise version of program. 

JET THRUST 

FJS 

Jet thrust, in newtons (lb), produced by jet 
lift engine when that option is used. Not 
included in cruise version. 


Multiple Cases 

The six data cards just described are necessary to determine the performance and 
geometry of a single gas generator. If a single case is to be run for each submission 
to the electronic computer, then both indexes N and NN on the first data card should be 
set equal to 1. However, the program is arranged so that several cases may be run per 
submission. Three methods are available for running multiple cases. In the following 
paragraphs, each method is described and illustrated with an example. 

The simplest method for running multiple cases is exercised when the input vari- 
ables to be changed are all among the six variables read from the sixth input data card. 
Then only one additional card is required for each case, and the desired values of the 
six variables are indicated on each card. The value of the index N on the first data card 
is set equal to 1, and the value of NN is set equal to the number of cases to be run. 
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For example, suppose it were required to investigate the effect of size on the 
characteristics of a gas generator. This could be done by changing the delivered gas 
flow WG. Suppose that the size range of interest could be covered by values of de- 
livered gas flow of 30, 35, and 40 kilograms per second. The order of data cards 
required to run these three cases would be as follows: 


Card 

Partial contents 

Card 

Partial contents 

1 

N = 1, NN = 3 

5 


2 


6 

WG = 30 

3 


6 

WG = 35 

4 


6 

WG = 40 


The second method is employed when all of the changes of input values are among 
the 14 variables appearing on the last two data cards. Then the first four data cards 
need be submitted only once, and they are followed by the appropriate combinations of 
the fifth and sixth cards. The same number of sixth cards must follow each fifth card. 
This number is the value given to the index NN. The index N takes on the value of the 
number of fifth cards to be submitted. 

For example, suppose that it were required to examine the effect of size (WG) on 
gas generators designed for a standard day (i. e. , TO = 288) and for a hot day (i. e. , 

TO = 305). The size effect would be studied by running the same three values of WG 
as in the first example, and for two different ambient temperatures, so that six cases 
would be required. The order of data cards would be as follows: 


Card 

Partial contents 

Card 

Partial contents 

1 

N = 2, NN = 3 

6 

WG = 35 

2 


6 

WG = 40 

3 


5 

TO = 305 

4 


6 

WG = 30 

5 

TO = 288 

6 

WG = 35 

6 

WG = 30 

6 

WG = 40 
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The third multiple -case method is used when the parameter to be changed appears 
on the first four data cards. Then a complete set of data cards has to be submitted for 
each value of this parameter. However, either of the first two methods may be com- 
bined with the third method. 

Suppose it was desired to investigate the effect of pressure loss in the inlet of the 
gas generator and that this effect was to be compared for two different ambient tem- 
perature designs, but for a single value of delivered gas flow. The effect of the inlet 
loss could be determined by changing the inlet recovery PI2P1 from a value of 0.99 
(assumed to be the value used for the examples under the first two methods) to 
PI2P1 = 0. 95. The required data input cards would be as follows: 


Card 

Partial contents 


Card 

Partial contents 

1 

N = 2, NN = 1 


1 

N = 2, NN-= 1 

2 



0 


3 

PI2P1 = 0.99 


3 

PI2P1 = 0.95 

4 



4 


5 

TO = 288 


5 

TO = 288 

6 

WG = 35 


6 

WG = 35 

5 

TO = 305 


5 

TO = 305 

6 

WG = 35 


6 

WG = 35 


Note that the data from all three examples could have been generated in a single 
submission by using the following set of data input cards: 


Card 

Partial contents 

Card 

Partial contents 

1 

N = 2, NN = 3 

6 

WG = 35 

2 


6 

WG = 40 

3 

PI2P1 = 0. 99 

1 

N = 2, NN = 1 

4 


O 






5 

TO = 288 

3 

PI2P1 = 0.95 

C 

V 

TTiri on 

VV KJ = 0\j 

A 

*1 


6 

WG = 35 

5 

TO = 288 

6 

WG = 40 

6 

WG = 35 

5 

TO = 305 

5 

TO = 305 

6 

WG = 30 

6 

WG = 35 
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COMPUTER PRINTOUT 


A typical sheet of computer printout for a single gas generator design is shown in 
figure 4. Figure 4(a) shows the printout in SI Units, and figure 4(b) gives the corre- 
sponding results in U. S. Customary Units. The first line on the page gives the name 
of the program and indicates whether the masses shown on the sheet correspond to lift 
or cruise engine technology. The title, data card 2, is printed next. 

All the inputs described in the section INPUT PARAMETERS follow. They are 
divided, somewhat arbitrarily, into primary and secondary inputs. 

If a turbine cooling air schedule has been used, the kind of schedule and the amount 
of cooling air are indicated on the next two lines following the inputs. 

The output from the program follows next. The first part of it is divided into two 
sections. The left-hand section contains parameters which describe the gas at the 
engine exhaust and parameters which pertain to the overall gas generator. The gas 
TEMPERATURE (TD) and PRESSURE (PD2) are those which prevail at the duct-nozzle 
exit. The units are K (°R) and kilonewtons per square meter (lb/sq ft), respectively. 
The specific POWER (GHP), or specific energy, is computed on the assumption that the 
gas expands isentropically to ambient pressure from the state at the end of the duct. 

The units are kilowatts per kilogram per second (hp/(lb/sec)). The SFC, specific fuel 
consumption (GSFC), is the ratio of the gas generator fuel flow to gas flow in kilograms 
of fuel per hour divided by kilograms of gas per second ((lb fuel/hr)/(lb gas/sec)). 

The DUCT DIAMETER (DUCTD) and STATIC PRES (PDSTAT) at the duct exit are 
given in meters (ft) and kilonewtons per square meter (lb/sq ft), respectively. Had the 
lift jet engine option been called for, these quantities would have applied to the nozzle 

v. 

exit. 

The remaining output variables in the left-hand section refer to the gas generator 
in its entirety. The FUEL FLOW (WF) is given in kilograms per hour (lb/hr). The 
JET THRUST (FJ) in newtons (lb) is computed on the assumption that the gas at the duct 
exit expands to ambient pressure through a nozzle with the velocity coefficient CFJ. 

The THRUST SFC (SFC) has the units of kilograms per hour per newton ((lb fuel/hr)/ 
lb thrust). The SPECIFIC THRUST (FJW1) is based on the gas generator inlet airflow, 
and is expressed in newtons per kilogram per second (lb/(lb/sec)). 

The right-hand section of the output lists the lengths and masses of the gas genera- 
tor components and exhaust duct and the percent of the total gas generator mass for each 
of the components. As indicated previously, because the size and mass of the inlet are 
so intimately related to the engine installation, no attempt was made to calculate a 
representative length or mass for this component. Lengths and masses of the gas 
generator and its components are in meters (ft) and kilograms (lb), respectively. The 
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length printed in parentheses under the total engine length is the length between the 
compressor inlet face and the turbine discharge. 

The characteristics of the turbine are printed on the next line of output. The data 
heading, FORTRAN symbol, and description of the parameter follows: 


NUMBER OF 
STAGES 

SNT 

Number of turbine stages. 

ALPHA1 = 
-BETA2 

ALPHAT 

Stator flow angle measured from the axis of rotation, 
in radians (deg). Because the velocity diagram is 
symmetrical, this angle is equal to the angle of re- 
lative flow leaving the rotor. 

BETA1 = 
-ALPHA2 

BETA 

Relative flow angle entering the rotor and absolute 
flow angle leaving the rotor, in radians (deg) . 

VX1 = VX2 
(= VO) 

VX 

Axial velocity through the turbine, in meters per 
second (ft/sec); also assumed to be the approach 
velocity to the first -stage stator. 

VU1 = 
-WU2 

VU1 

Stator tangential velocity component and tangential 
component of relative velocity out of the rotor, in 
meters per second (ft/sec). 

VU2 = 
-WU1 

VU2 

Tangential component of absolute velocity leaving the 
rotor and tangential component of the relative veloc- 
ity entering the rotor, in meters per second (ft/sec). 

STAGE 

LAMDA 

STLAMT 

Stage speed-work parameter, U^/(AH). 

MEAN 

SPEED 

UTM 

Mean blade speed U m , in meters per second 
(ft/sec). 

ABS INLET 
MACH NO. 

AMT1 

Mach number of the absolute velocity into the first- 
stage rotor. 

REL. INLET 
MACH NO. 

AMT1W 

Mach number of relative velocity at the first -stage 
rotor inlet. 
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FLOW 

FLOWT 

Ratio of through-flow velocity to mean blade speed. 

COEFF. 



ABS OUTLET 

AMT2 

Mach number of absolute velocity at discharge of 

MACH NO. 


last -stage rotor. 


The remaining output format is virtually self-explanatory. It displays the thermo- 
dynamic properties of the working fluid throughout the gas generator, as well as the 
principal dimensions of the components. 

The data are arranged in tabular form with the name of the component appearing in 
the first column. Parameters such as length (m (ft)), pressure ratio, change in enthal- 
py (kj/kg (Btu/lb)), fuel-air ratio of the working fluid, and total efficiency, all of which 
are applicable to the component as a whole, are printed on the same line as the name of 
the component. Values of parameters which are different for the inlet and outlet of the 
component are printed on lines above and below the component name line, respectively. 
These parameters include working fluid mass flow rate (kg/sec (lb/sec)), total temper- 
ature (K (°R)), total pressure normalized by standard sea-level pressure, axial Mach 
number, axial velocity of the working fluid (m/sec (ft/sec)), hub-tip ratio, and the 
corresponding tip and hub diameters (m (ft)). The hub diameter is printed within paren- 
theses under the corresponding tip diameter. 

The two sets of values for mass flow and temperature appearing at the inlet of the 
turbine are for the stator and rotor inlet, respectively. The values reflect the contri- 
bution of the stator cooling air, which is 0. 5 of the total turbine cooling air. The re- 
maining half of the cooling air is added downstream of the turbine. 

In the column headed PRESS RATIO, the value of the ratio is determined by dividing 
the outlet pressure by the inlet pressure. For the turbine, the reciprocal of the pres- 
sure ratio, as defined previously, is also printed within parentheses. 

The last number on the table, under TIP DIA. , is the diameter in meters (ft) of the 
exit of the exhaust duct. If the jet lift engine option had been used this last number 
would have been the exhaust diameter of the jet nozzle. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 17, 1972, 

501-14. 
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PROGRAM 4L LIFT ENGINE ROTATING LUnpuNtNii 
ONE SPOOL GAS GENERATOR/LIFT ENGINE 

DIMENSION TITLE172) 

COMMON / JANI/ AJ,PIE,G,ZZ f PSTO, TSTO, PS I R , RA, DELO, SOTH I 
DIMENSION UNIA(4), UNI0I6) 

DATA UNIA(l) , UNI A ( 2 ) »UNIA(3) ,UNIA(4)/6H CUS , 6HT0MARY , 6H , 6H 

1 SI/ 

DATA (UNI0(I),I=1,6)/6H G.6HAS GEN, 6HER ATOR , 6H , 6H LIFT , 

16HENGINE/ 

READ (5,58) JCGEOM, K I NO, N, NN 
READ (5,60) (TITLE)! ), 1-1,12) 

READ (5,56) ALPH AT, AM02, PC A , OHDTC l , AKCT , HF , P I2P 1 , AMC 1 
REAU (5,56) VC2VC1,ETAB, AKUT, UNITS, BUSER.CFJ 
IF (UNITS. EQ.O. ) GO TO 4 
ALPHAT=ALPHAT/. 01745 
HF=HF/ 2. 326 
DO 54 J=l,N 

READ (5,56) ETAC , DTMDC 1 , UT I PCC, PO , TO, PB2PB 1 , SNC , PD2PD l 
IF (UNITS. EQ.O. ) GO TO 6 
D I V=. 3048 

utipc:=utipcc/div 

PO=PO/. 04788 
T0= TO/ . 5556 
DO 52 K=1,NN 

READ (5,56) WG , PC2PC 1 , TT l , SNT , F JS , OPT I ON 
IF ( WG.EQ.O. 0} WG=50. 

IF (UNITS. EQ.O. ) GO TO 8 
WG=WG/ .4536 
TTl=TTl/.5556 
TCI =T0 

HC1=P0LY(TC1,1 ) 

OELO=PO/PSTD 
DELCI=P0*PI2P1/PSTD 
SQTHl=SQRT( TCI/519. ) 

CPC l=POLY ( TC I » 3 ) 

GAMCl=l./(l.-(RA/(AJ*CPCl) ) ) 

TEM 1=G AMC 1—1. 

TEM2=TEMl/2. 

TEM4=1.+(TEM2*AMC1*AMC1) 

TEM5=SQRT (GAMC1*G/RA ) 

TEM6=(GAMC1+1.)/(2.*TEM1 ) 

SWC1=TEM5/ZZ*AMCI/(TEM4**TEM6)*( 1 .-DHDTC1**2 ) 
VC1=AMC1*SQRT(GAMCI*G*RA*TCI/TEM4) 

CALL COMPRS ( TC 1 , PC2 PC 1 , DELC 1 , ET AC, HC1 , SWC l , TC2 , DELC2, DELHC , HC2 ) 

PC1=DELC1*PSTD 
PC2=DELC2*PST0 
CPC2=P0LY(TC2,3) 

G AMC2= 1 • / ( 1 .- ( R A/ ( AJ*CPC 2) ) ) 

VC 2= VC 1*VC2VC l 

AMC2=SQRT(VC2*VC2/(GAMC2*G*RA*TC2-(GAMC2-1. )/2.*VC2*VC2) ) 
RH0C2=DELC2*PST0/(G*RA*TC2)*( 1 .-VC2**2/ ( 2 .* A J*G*CPC2*TC2 ) )**( 1./ (G 
IAMC2-1. ) ) 

TEM9=G*RH0C2*VC2 
TC00L=TC2 

HCOOL=POLY(TCOOL,1) 

W1AC1=SWC1*DELC1/S0TH1 

IF JCGEOM = 1 CONSTANT HUB DIAMETER 

IF ( JCGEOM. EQ.l) GO TO 10 

IF JCGEOM = 2 CONSTANT MEAN DIAMETER 

IF ( JCGEOM. EQ. 2) GO TO 12 
C 

C IF JCGEOM = 3 CONSTANT TIP DIAMETER 


A 

1 


A 

2 


A 

3 


A 

4 


A 

5 


A " 

6 


A 

7 


A 

8 


A 

9 


A 

10 


A 

11 

1 

A 

12 

6 
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13 

13 

A 

14 

14 
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15 


A 

16 


A 

17 


A 

18 


A 

19 

21 

A 

20 


A 

21 


A 

22 


A 

23 


A 

24 


A 

25 


A 

26 

27 

A 

27 


A 

28 


A 

29 


A 

30 


A 

31 


A 

32 

34 

A 

33 


A 

34 


A 

35 

35 

A 

36 

36 

A 

37 


A 

38 


A 

39 


A 

40 


A 

41 

37 

A 

42 


A 

43 

38 

A 

44 


A 

45 

39 

A 

46 


A 

47 

40 

A 

48 


A 

49 


A 

50 

42 
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51 


A 

52 


A 

53 

43 

A 

54 


A 

55 

44 

A 

56 
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A 

58 
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Figure 3. - FORTRAN listings of gas generator/lift engine computer program. 
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c 


10 


12 


14 

16 


C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 


18 

20 

22 


IF ( JCGEOM.EQ. 3) GO TC 14 

0H0TC2=1. /SQRTI 1.+W1ACI/ < TEM9*DHDTC1**2 ) ) 

DC2DCI=DHDTC1/DHDTC2 
GO 10 16 

TEM L 0= ( (1.+DH0TCD/2. )**2 
TQ=WlA:i/ <4.*TEM9*TEM10) 

DHDTC2=(1 .-TQ)/( l.+TQ) 

0C2DC l = ( l . +DHDTC 1 ) / ( 1.+DHDTC2) 

GO TO 16 

OHO TC2 = SORT < 1.-W1AC1/TEM9) 

0C2DC1 = 1. 

DEL82=DELC2*PB2PB1 
DELT 1=DELB2 
HT1=P0LY(TT1,1) 

PSIHT1 =POLY( TT1 ,4) 

HO=HF+ 1254. 

F8= ( HT 1-HC2 ) / ( HO- (HC2/ETAB)-PSIHTl) 

FAR=F3/(£TAB*(1.-FB/ETAB>) 

HTIG ENTHALPY OUT OF COMBUSTOR 
CHRGB CHARGEABLE BLEED 
HTIGT EFFECTIVE ENTHALPY AT ROTOR INLET 

HN EFFECTIVE ENTHALPY AT TURBINE EXHAUST AFTER COOLING AIR IS MIXED 
DELHT TURBINE ENTHALPY DROP DUE TO WORK REQUIRED BY COMPRESSOR 
AND PUMPING OF COOLInu AIR 
BUSER = BLEED FLOW! LBS/SEC ) /W 1 
IF KIND= 0 PROGRAMMER FURNISHES COOLING BLEED 
IF KIND = 50 LIFT ENGINE 
IF KIND = 100 CRUISE ENGINE 

I NOE X=K I ND/50 
IF (INDEX-1) 18,20,22 

pca=p;a 

GO TO 24 

PCA=.000110*TT1— .242 
GO TO 24 

PC A =. 0001 50* TT 1—. 297 


A 70 
A 71 
A 72 
A 73 
, 74 
Ik 75 
A 76 
A 77 
A 78 
A 79 
A 80 
A 81 
A 82 
A 83 
A 84 
A 85 
A 86 
A 87 
A 88 
A 89 
A 90 
A 91 
A 92 
A 93 
A 94 
A 95 
A 96 
A 97 
A 98 
A 99 
A 100 
A 101 
A 102 
A 103 
A 104 
A 105 
A 106 
A 107 


24 


26 

C 


C 

C 

C 


28 


30 


CHRGB= • 50 

BETCOL=PC A* ( 1. -BUSER) 

BET TOT = BE TCOL+ BUSER 
TEMll=< 1. +FAR)*( l.-BETTOT) 

HTIG=HT1+FB*PSIHT1 
DELHB=HTI G— HC2 
UTM=0TM0CI*UTIPCC*SQTH1 
DHPUMP=UTM*UTM/AJ/G 

OELHT= ( DE LHC+DHPUMP*CHRGB*BETCOL ) / ( TEM1 l + ( l .-CHRGB ) *BETCOL ) 

HCOLP=HCOOL*DHPUMP 

W1=WG/ (TEM1H-BETCQL) 


WT=WG 

WB1=W1*< l.-BETTOT) 

WB2=WB1*( l.+FAR) 

AC1 = WI*SQTH1/(DELC1*SWC1 ) 

DC 1 = SQRT ( 4 . * AC1 /P I E ) 

DTM=DC1*DTMDC1 

ROTOR INLET CONDITIONS 

HT I GT= ( HT I G*TEM1 1+ ( 1 . -CHRGB ) *8ETC0L*HC00L ) / ( TEM1 1 ♦( 1 .-CHRGB ) *BETCO 
IL) 

FARR=FAR/(1 .+ BET COL* ( l . -CHRGB ) / ( l . -BETTOT ) ) 

FR=ETAB*FARR/(1.+FARR> 

T TR=TT 1 

DELTAT=(HTIGT-P0LY(TTR,l)-FR*P0LY(TTR,4) )/(POLY(TTR,-l)+FR*POLY(TT 
IR ,- 4 ) ) 

TTR=TTR+OELTAT 

IF (ABS(OELTAT).LT.. 1 ) GO TO 30 
GO TO 28 

WR=W1*(TEM11+BETC0L*( l. -CHRGB ) ) 


A 108 
A 109 
A 110 
A 111 
A 112 
A 113 
A 114 
A 115 
A 116 
A 117 
A 118 
A 119 
A 120 
A 121 
A 122 
A 123 
A 124 
A 125 
A 126 
A 127 
A 128 
A 129 
A 130 
A 131 
A 132 
A 133 
A 134 
A 135 
A 136 
A 137 
A 138 
A 139 


56 


61 


63 

64 


74 


76 77 78 79 


Figure 3. - Continued. 
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c 


A 

140 




CALL TURBIN (TTR,DELT1,DELHT, AKCT *AKUT,WR, ALPHAT »SNT,FR, l.,DCl,UTM 

A 

141 




1 ,0TMDC1 ,TT2,PT2PT1,DELT2,AMT2X,0HDTT2, AMBBT, OT 2DC 1 , ET ABT, VX, VU1, VU 

A 

142 




22,BETA,0H0TTI,DT10C1,AMTIX,AMTIW,AMT1, FLOKT) 

A 

143 



c 


A 

144 

85 



PT1PT2=1./PT2PT1 

A 

145 




stlamt =SNT*AM8BT 

A 

146 




BETAR=BETA*.01745 

A 

147 




AMT2=AMT2X/COS(8£TAR ) 

A 

148 

87 



HN=( ( HTI GT-DELHT ) *(T EMl l + ( 1 . -CHRGB ) *BETCOL ) +CHRGB* BETCOL*HCOLP > / ( T 

A 

149 




1EM1H-8ETCOL) 

A 

150 




FAPRlM=FAR/< l.+BETCOL/I 1 .-BETTOT) ) 

A 

151 




fprim£=etab*faprim/< l.+FAPRIMI 

A 

152 




TD=TT2 

A 

153 



32 

0£L T = t HN— (P0LY(TD,l)+FPRIME*P0LY(TD,4) ) > / < POL Y ( T D, - L ) *FPR I ME*POL Y { 

A 

154 




l TO » -4) ) 

A 

155 

89 

90 


TD=TD*-D£LT 

A 

156 




IF ( aBS(DELT) .LT..X) GO TO 34 

A 

157 




GO TO 32 

A 

158 



C 


A 

159 



C 

BEGIN DUCT CALCULATIONS 

A 

160 



C 


A 

161 



34 

CPD=P0LY( TO,3)+FPRIME*POLY(TD,6) 

A 

162 

98 

99 


rd=ra*fprime*psir 

A 

163 




GAMMD^l./I l.-RD/(AJ*CPD) ) 

A 

164 




WDW1=HG/Wl 

A 

165 




P01=D£LT2*PSTD 

A 

166 



C 


A 

167 




CALL DUCT 2 ( PD1, P02PD1. AM02,TD,WDW1 , SWC L , OELC 1 , GAMMO, RD, PDSTAT , AD2 

A 

168 




LAI ) 

A 

169 



C 


A 

170 

100 



PD2=P?1*PD2PD1 

A 

171 




VJ=SQRT(2. *GAMMD/ { GAMMO— L . ) #G*RD*TO*I l • - ( P0/P02 ) ** ( ( GAMMD— 1 • ) /GAMM 

A 

172 




ID) ) )*CFJ 

A 

173 

102 

103 


DEL02*=PD2/PST0 

A 

174 




DVJCTD^SQRT 10C1*DC1*AD2 A1 » 

A 

175 



c 


A 

176 

104 



VD2=S3RT(GAMMD*G*R0*T0*AMD2*AMD2/< 1 .0+1 GAMMD-l . 0 ) *AMD2*AMD2/2.0 ) ) 

A 

177 

105 



FJ=WOrfl*WL*VJ/G 

A 

178 




IF (OPTION. LT. 1.5) GO TO 36 

A 

179 




OF J=F J -F JS 

A 

180 




FJTESr=OFJ/FJS 

A 

181 




IF (A)S(FJTEST).LT..001) GO TO 36 

A 

182 




WG=FJS*G/VJ 

A 

183 




GO TO 26 

A 

184 



36 

AMN=SflRT( VJ*VJ/( GAMMD*G*RO*TO-( GAMMD- 1.0 )/2.*VJ*VJ ) ) 

A 

185 

114 



WF=WlPFAR*( I. -BETTOT )*3600. 

A 

186 




SFC=WP/FJ 

A 

187 



C 


A 

188 



C 

3 AS ANO THRUST SFC'S GAS POKER 

A 

189 



c 


A 

190 




gsfc=hf/wg 

A 

191 




SFCNTH=WF/FJ 

A 

192 




FJW1=)IDW1*VJ/G 

A 

193 




PHI X=POLY (TD»2 ) + FPRI ME*POLY( TO, 5 ) 

A 

194 

115 

116 


PHI XI* PHI X-RD/A J*ALOG ( PD2/ PO ) 

A 

195 

117 



TXI -TO 

A 

196 



38 

DELTX*TXI*(PHIXI-P0LY(TXl,2)-FPRIME*P0LY(TXI,5) ) / 1 POL Y ( TX 1 , 3 )+FPR I 

A 

197 




1ME*P0UY(TXI ,6) ) 

A 

198 

119 

120 


txi=txi*deltx 

A 

199 




IF (D6LTX.GT..1) GO TO 38 

A 

200 




HXI=P0LY(TXI ,1 )+FPRIME*POLY< TXI, 4) 

A 

201 

126 

127 


HD=P0LY(TD,1)+FPRIME*P0LY(T0,4) 

A 

202 

128 

129 


GHP=AJ*(H0-HXI 1/550. 

A 

203 



C 


A 

204 



C 

COMPONENT 01 AMETERS 

A 

205 



C 


A 

206 




DC 1 T=3C 1 

A 

207 




DC1H=0C1*DHDTC1 

A 

208 




DC2 T=‘3C 1 T *0C20C 1 

A 

209 




0C2H=DC2T*DHDTC2 

A 

210 




Figure 3. - Continued. 
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ooo noo ooo o o o o o o noo ooo non noo o o o o o o o 


DT1T=3C1T*DT10C1 A 211 

DT1H=0UT*DH0TTI A 212 

DT2T=3C1T*DT2DC1 A 2 1 B 

DT2H=3T2T*0HDTT2 A 214 

A 215 
A 216 

LENGTHS - LIFT ENGINE ROTATING COMPONENTS A 217 

A 218 
A 219 

COMPRESSOR A 220 

A 221 

CLENDM*0.2+< . 234-. 2T8*0HDTC1 >*SNC A 222 

C0MPL=CLENDM*(DClH»DClT)/2. A 223 

A 224 

COMBUSTOR A 225 

A 226 

WBIN=Wl*(l.— BETTOT) A 227 

CORWBl =WBIN*SQRT (TC2/TST0I/DELC2 A 228 130 

EL0H=2 . 0 A 229 

VREF=80. A 230 

DMAV=< QC2T*DC2H+DTlT+0TlH)/4. A 231 

8URNL = 2 . 2*CQRW81 *EL0H*SQRT I TC2 ) / ( DMAV*VREF ) / 1 2 • A 232 

A 233 

TURBINE A 234 

A 235 131 

TSAR=5.45-5.97*DHu.Tl A 236 

TRAR=5.1-5.5*0H0TT1 A 237 

IF (TRAR.GT.6.) TRAR=6. A 238 

BLDH=3UT*< 1.-0HDTT1 )/2. A 239 

SAXCHD=BLDH/TSAR A 240 

RAXCH0=BLOH/TRAR A 241 

TCLR=. A0ORAXCHD A 242 

TUR8L = SNT * ( SAXCHO+RAXCHO+2. *TCLR ) +TCLR+SAXCHD A 243 

A 244 

DUCT A 245 

A 246 

DUCTL=.75*DT2T A 247 

IF I0PTI0N.GT.1.5) DUCTL=.25*DT2T A 248 

A 249 

ENGINE A 250 

A 251 

ENGL=C0MPL+BURNL+TURBL A 252 

ENGLT=ENGL+DUCTL A 253 

EL0D=ENGL/DC1 A 254 

EL00T=ENGLT/0C1 A 255 

A 256 

WEIGHTS - LIFT ENGINE ROTATING COMPONENTS A 257 

A 258 
A 259 

COMPRESSOR A 260 

A 261 

ELOORF = . 2+. 08 1*SNC A 262 

Q=.5+.54CLENDM/EL0DRF A 263 

C0MPW= 5.0*( .25*( DC1 T ♦ DC1 H+DC2T +DC2H ) ) **2 . 2 *SNC**1 . 2 ( UT I PCC*SQTH A 264 

11/1103. )**. 3 A 265 

A 266 

COMBUSTOR A 267 

A 268 136 137 

DMBURN= ( 0C2T ♦DC2H+DT IT +0T 1H ) / 4. A 269 

8URNW= 32 . *DMBURN*DMBURN*SQRT ( ELOH ) A 270 

A 271 

TURBINE A 272 

A 273 139 

DMTURB = (DTI T+DT1H+DT 2T+DT2H ) /4. A 274 

TURB W=. 264SNT 4DMTURB **2 .5*UTM**0.6 A 275 

A 2 76 

DUCT A 277 

A 278 140 141 

0UCTW=4.7*DT2T*0T2T A 279 

IF (OPTION. GT. 1.5) DUCTW=1 . 1 *0T2T*0T2T A 280 
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ooo o ooo ooo 


A 281 

ACCESSORY. STRUCTURE, AND TOTAL A 282 

A 283 

ACCW=» 002* ( F J+l . 35*WF ) A 284 

8ENGW=C0MPW+BURNW*TURBW*ACCW A 285 

STRW=. 10*BENGW A 286 

ENGW=8ENGW+DUCTW+STRW A 287 

PCTH=100./ENGH A 288 

C0MPWX=C0MPW*PCTW A 289 

8URNWX=8URNW*PCTW A 290 

TURBWX=TURBW*PCTW A 291 

DUCTWX=OUCTW*PCTW A 292 

ACCWX=ACCW*PCTW A 293 

STRWX=STRW*PCTW A 294 

AJCGE0- JCGEOM A 295 

AKIN0=KIND A 296 

IF (UNITS. EQ. 0. ) GO TO 40 A 297 

A 298 

WRITE INPUTS IN SI UNITS A 299 

A 300 

WG=.4536*WG A 301 

TT1=.5556*TT1 A 302 

P0= . 047884P0 A 303 

TO=.5556*TO A 304 

A 305 

UTIPCC=DIV*UTIPCC A 306 

HF=2.326*HF A 307 

ALPHAT =.01745*ALPHAT A 308 

A 309 

WRITE OUTPUTS IN SI UNITS A 310 

A 311 

TD=.5556*TD A 312 

P02=.04788*PD2 A 313 

GHP=l.644*GHP A 314 

GSFC=SSFC A 315 

DUCTO=DUCTD*« 3048 A 316 

POSTAT=PDSTAT*. 04788 A 317 

WF=.4536*WF A 318 

FJ=4.448*FJ A 319 

SFC=.1020*SFC A 320 

FJW1=9.806*FJWI A 321 

C A 322 

WTC0N= .4536 A 323 

COMPW=WTCON*COMPW A 324 

BURNW=WTCON*BURNW A 325 

TURBW= WTC 0N*TURBW A 326 

ACCW=WTCON*ACCW A 327 

STRW=WTCCN*STRW A 328 

ENGW=WTCGN*ENGW A 329 

C A 330 

BETA=. 01745*BETA A 331 

VX=DIV*VX A 332 

VU1=DI V*VU1 A 333 

VU2=DIV*VU2 A 334 

UTM=0IV*"TH A 335 

c A 336 

ELC0N=,3048 A 337 

C0MPL=ELC0N*CQMPL A 338 

8URNL=ELC0N*8URNL A 339 

TURBL=ELC0N*TURBL A 340 

0UCTL=EUC0N*0UCTL A 341 

ENGL=ELCON*ENGL A 342 

ENGLT=ELC0N*ENGLT A 343 

C A 344 

HCON=2.326 A 345 

OELHC=HCON*DELHC A 346 

OeLHB=HCON*U£LHB A 347 

OELHT=HCON*DELHT A 348 

C A 349 

W1=WTC0N*W1 A 350 

W81=WTC0N*WB1 A 351 
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WB2=WTC0N*WB2 

A 

352 



WR=WTC 

ON*WR 

A 

353 



WT=WTC 

ON*WT 

A 

354 


c 



A 

355 



TCON=. 

5556 

A 

356 



TC1=TC 

ON* TCI 

A 

357 



tc2=t: 

0N*TC2 

A 

358 



ttr=t: 

ON*TTR 

A 

359 


c 



A 

360 



VCl =EL C0N*VC1 

A 

361 



VC2 =EL 

C0N*VC2 

A 

362 



VD2=ELC0N*VD2 

A 

363 


c 



A 

364 



DC1T=ELCGN*DCIT 

A 

365 



CC1H=ELC0N*DC1H 

A 

366 



DC2T=ELCCN*0C2T 

A 

367 



0C2H=E 

LCCN*DC2H 

A 

368 



DT1T=ELCGN*DT1T 

A 

369 



DT1H=ELC0N*DT1H 

A 

370 



DT2T=ELCGN*DT2T 

A 

371 



0T2H=ELCGN*DT2H 

A 

372 


c 



A 

373 


40 

WRITE 

(6,62) 

A 

374 

147 


WRITE 

(6,60) (TITLE! I ), 1=1, 12) 

A 

375 

148 


WRITE 

(6,110) 

A 

376 

153 


WRITE 

(6,72) WG,PO,UTIPCC,ALPHAT 

A 

377 

154 


WRITE 

(6,74) TT1,T0,SNC,AKCT 

A 

378 

155 


WRITE 

(6,76) PC2PC1,ETAC,0TM0C1 

A 

379 

156 


WRITE 

(6,78) SNT, AKI ND, A JCGEO, AKUT 

A 

380 

157 


IU=1 


A 

381 



IF (UNITS. NE. 0.0) I U= 3 

A 

382 



WRITE 

(6,80) UNIA( IU) ,UNIA( I U+l ) , P I 2P 1, VC2VC 1 

A 

383 

160 


10=1 


A 

384 



IF (OPTION. GT. 1.5) 10=4 

A 

385 



IOE=IO+2 

A 

386 



WRITE 

(6,82) (UNI0(IP),IP=I0,I0E), AMCl.HF.CFJ 

A 

387 

166 


WRI TE 

(6,84) F JS , OHDTC 1,ETAB,P02PD1 

A 

388 

171 


WRI TE 

(6,86) BUSER.PB2PBI, AMD2 

A 

389 

172 


IF (KIND.EO.O) GO TO 42 

A 

390 



IF (KINO. EQ. 50) GO TO 44 

A 

391 



IF (KINO. EO. 100) GO TO 46 

A 

392 


42 

WRITE 

(6,64) 

A 

393 

182 


WRI TE 

(6,70) PC A 

A 

394 

183 


GO TO 

48 

A 

395 


44 

WRITE 

(6,66) 

A 

396 

185 


WRI TE 

(6,70) PC A 

A 

397 

186 


• GO TO 

48 

A 

398 


46 

WRITE 

(6,68) 

A 

399 

188 


’ WRITE 

(6,70) PC A 

A 

400 

189 

48 

WR I TE 

(6,112) 

A 

401 

190 


WRI TE 

(6,88) 

A 

402 

191 


WRITE 

(6,90) TO 

A 

403 

192 


WRITE 

(6,92) P02,C0MPL,CGMPW,C0MPWX 

A 

404 

193 


WRI TE 

(6,94) GHP , BURNL , BURNW , BURNWX 

A 

405 

194 


WRITE 

(6,96) GSFC,TUR8l,TURBW,TURBWX 

A 

406 

195 


WRI TE 

(6,98) DUCTO.OUCTL, OUC TW , DUCT WX 

A 

407 

196 


WRI TE 

(6,100) POST AT , ACCW, ACCWX 

A 

408 

197 


WRITE 

(6,102) WF , STR W , STRWX 

A 

409 

198 


WRITE 

(6,104) FJ 

A 

410 

199 


WRITE 

(6,106) SFC.ENGLT, ENGW 

A 

411 

200 


WRITE 

(6,108) FJW1.ENGL 

A 

412 

201 


WRI TE 

(6,118) 

A 

413 

202 


WRITE 

(6,116) 

A 

414 

2 03 


WRITE 

(6,114) SNT,ALPHAT,BETA,VX,VU1, VU2, STL AMT , UTM, AMT1 , AMT1 W, F LO 

A 

415 



1WT, AMT2 

A 

416 

204 


WRI TE 

(6,120) 

A 

417 

205 


WRITE 

16,136) 

A 

418 

206 


WRITE 

(6,134) W1,TC1 ,0ELC1,AMCI,VC1, DHDT Cl , 0C1T,0C1H 

A 

419 

207 


WRITE 

(6,124) PC2PC1 .DELHC.ETAC 

A 

420 

208 


WRI TE 

(6,126) WB1,TC2,DELC2,AMC2,VC2, DHDTC2, DC2T,DC2H 

A 

421 

209 
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50 

52 

54 

C 

56 

58 

60 

62 

64 

66 

68 

70 

72 


74 


76 

78 


80 

82 

84 


86 

88 

90 

92 

94 

96 

98 

100 

102 

104 

106 

108 

110 

112 

114 

116 


WRITE (6,122) PB2PB1 , DELHB , F AR, ET AB 

WRITE (6,132) WB2,TT1,0ELT1, AMT1X, VX, DHDTT 1 , OT 1 T , WR , TTR, OT 1H 
WRITE (6,128) PT2PT1 , DELHT , F ARR, ET ABT 

WRITE (6,130) PT1PT2,HT«TD,DELT2,AMT2X,VX,DHDTT2»DT2T »DT 2H 
AMD2S=AMD2 

IF (OPTION. LE. 1.5) GO TO 50 

AMD2=AMN 

VD2=VJ 

WRITE (6,138) PD2PD1 ,FAPRIM,WG»TD, DEL 02 , AHD2 , V02, DUCT 0 

AMD2=AM02S 

CONTINUE 

CONTINUE 

GO TO 2 


, F8 . 3, 24H 
• F8.3.24H 


AMBIENT PRESSUR 
TURB NOZZLE ANGL 


,F8.3,24H 
, F8 . 3, 24H 


AMBIENT TEMP 
TURB LOSS COEF 


FORMAT (8F10.0) 

FORMAT (415) 

FORMAT (12A6) 

FORMAT (1H1,45X,5 IHONE SPOOL GAS GENERATOR/LIFT ENGINE-LIFT COMPON 
1 ENT S) 

FORMAT (IX, 16H COOLING BLEED) 

FORMAT (IX.11HLIFT ENGINE) 

FORMAT (IX, 13HCRUISE ENGINE) 

FORMAT (IX, 9HPC A =,E12.5) 

FORMAT ( 1H0.24H GAS FLOW 
IE ,F8.3,24H CORR TIP SPEED 

2E ,F8.3) 

FORMAT ( 25H TURBINE INLET TEMP 

1 , F8 . 3, 24H NO. COMP. STAGES 

2 , F8 . 3 ) 

FORMAT ( 25H COMP PRES RATIO , F8 . 3, 32X , 24H COMP EFFICIENCY 

l , F8. 3 ,24H FLARE RATIO! 0TM/0C1 ) ,F8.3> 

FORMAT (25H TURBINE STAGES ,F8.3,24H ENGINE APPLICATION 

1 ,F8.3,24H COMP FLOW PATH ,F8.3,24H TURB STRAIGHT VANES 

2 ,F8.3) 

FORMAT (21H UNITS ,2A6,24H INLET PRES RATIO ,F 

18.3.24H AXIAL VELOCITY RATIO ,F8.3) 

FORMAT (15H OPTION .3A6.24H INLET AXIAL MACH NO ,F8.3,23 

IH FUEL HEATING VALUE ,F9.3,24H NOZZLE COEF ,F8.3) 

FORMAT (25H JET THRUST ,F8.3,24H INLET HUB-TIP RATI 

10 , F3 . 3 , 24H COMBUSTOR EFFICIENCY ,F8.3,24H DUCT PRES RATIO 

2 , F8 . 3 ) 

FORMAT (33X.24H USER BLEED ,F9.4,23H COMBUSTOR PRES RA 

1TI0 , F 8. 3 , 24H DUCT EXIT MACH NO ,F8.3) 

FORMAT (1H0.13X, 14HGAS PROPERT I ES , 57X , 6HLENGT H, 7X , 4HMASS, 6X , 7HPE RC 
1 ENT ) 

FORMAT (20X, 1 lHTEMPERATURE.Fl 1.1 ) 

FORMAT (20X, 11HPRESSURE , F 1 1 . 1 , 25X , 13HC0MPRESS0R 


,F11.1,25X, 13HC0MBUST0R 
, Fll.l,25X,l 3HTURB INE 


,Fl 1.3, 2FU.1 
,Fll.3,2Fll.l 
, F 1,1. 3, 2F 11. 1 


,F13.3,23X, l 3HDUCT /NOZZLE ,F11.3, 


, 11X.2F 


, 11X.2F 


1) 

FORMAT ( 20X , 1 IHPOWER 
1) 

FORMAT (20X.11HS F C 
1 ) 

FORMAT ( 14X, 17HDUCT OIAMETER 
12F11.1 ) 

FORMAT (14X.17H STATIC PRES , F 1 l . 1 , 2 5X , 1 3HACCESS0R I ES 

111.1) 

FORMAT ( 1 4X , 1 7HFUEL FLOW 
111.1) 

FORMAT (14X.17HJET THRUST 

FORMAT ( 14X , 17HTHRUS T SFC 
1 ,6X,5H100.0) 

FORMAT (14X, 17HSPECI FIC THRUST 

FORMAT (IH0,9X»1 4HPR I M AR Y I N PUTS , 5 IX , 16HSEC0ND ARY INPUTS) 

FORMAT ( 1H0,57X,6H0UTPUT ) 

FORMAT ( 1X,7HTUR8INE,9X,F3.0,5X,F7.2,2X,F7.2,2X,F7. I , IX , F7 . 1 , FB. 1 , 
11X,F7.4,1X,F7.1,2X,F7.4,5X,F7.4,3X,F6.3,4X,F6.3//) 

FORMAT (1H0, 14X, l 09H NUMBER OF ALPHA1= 8ETA1= VX1=VX2 VU1= 
IVU2= STAGE MEAN ABS. INLET REL. INLET FLOW ABS OUTLET/16 
2X, 109HSTAGES -BETA2 -ALPHA2 !=V0) -WU2 -WU1 L AMD A 

3SPEED MACH NO. MACH NO. COEFF. MACH NO. ) 


,Fll.l,25X,13HSTRUCTURE 

, F 1 1 . 1 ) 

,F13.3,27X,9HT0TAL , F 1 1 . 3, F ! 1 . 1 

,F11.1,43X,1H( »F5.3, 1H) ) 


A 422 
A 423 
A 424 
A 425 
A 426 
A 427 
A 42 8 
A 429 
A 430 
A 431 
A 432 
A 433 
A 434 
A 435 
A 436 

A 438 
A 439 
A 440 
A 441 
A 442 
A 443 
A 444 
A 445 
A 446 
A 447 
A 448 
A 449 
A 450 
A 451 
A 452 
A 453 
A 454 
A 455 
A 456 
A 457 
A 458 
A 459 
A 460 
A 461 
A 462 
A 463 
A 464 
A 465 
A 466 
A 467 
A 468 
A 469 
A 470 
A 471 
A 472 
A 473 
A 474 
A 475 
A 476 
A 477 
A 478 
A 479 
A 480 
A 481 
A 482 
A 483 
A 484 
A 485 
A 486 
A 487 
A 488 
A 489 
A 490 
A 491 


210 

211 

212 

213 


217 
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118 FORMAT (1HO) A 492 

120 FORMAT ( 1H0 , 22X, 5HPRE SS , 14X , 88HFUEL- A I R TOTAL MASS TOTA A 493 

1L TOTAL AX. MACH AXIAL HUB-TIP TIP D I A. / 12X, 1 19H A 494 

2 RATIO DELTA H RATIO EFF. FLOW TEMP A 495 

3 PRESS NUMBER VELOCITY RATIO (HUB DIA.)//) A 496 

122 FORMAT (10H COMBUSTOR , 8X , F 10 .4 , F9 . 2 , F 10 . 5, F 1 0. 5 ) A 49? 

124 FORMAT (10H COMPRESS , 8X ,F 10 . 4, F9 . 2 , 12X , F8 .5 ) A 498 

126 FORMAT ( 62 X, F7 . 2 , 3X, F7. 1 , F 10 .4, F9 . 4, F 1 0 . 1, F 10. 5 , FI 1 .4/ 120X , 2H (,F7 A 499 

1.4, 2Hj ) A 500 

128 FORMAT (10H TURBINE , 8X ,F 10 . 4, F9 . 2, F 10.5 , F 10. 5) A 501 

130 FORMAT (16X.5H (,F7.4,2H) , 32X , F7 . 2, 3X, F7 . 1 , F 10 .4, F9. 4, F 1 0. 1, F 1 A 502 

10.5,Fll.4/120X,2H (,F7.4,2H )) A 503 

132 FORMAT (62X , F7 . 2 , 3X, F7. 1 , F 10 .4, F9 . 4, F 10 . 1 , F 10 . 5 , F 1 1 . 4/62X, F 7 . 2 ,3X , A 504 

1F7.1»41X,2H (.F7.4.2H) ) A 505 

134 FORMAT 1 62X , F7. 2 , 3X, F7 . 1 ,F 10 .4, F9 . 4, F 10. 1 , FI 0 . 5, F 11 .4/ 120X , 2H (,F7 A 506 

1.4.2H) ) A 507 

136 FORMAT (12H INLET > A 508 

138 FORMAT (12H DUC T/NOZ ZLE , 6X , F 10 . 4, 9X, F 10. 5, 1 5X , F 7 . 2 , 3X, F7 . 1 /79X, F 1 0 A 509 

1.4,F9.4,Fl0.l,10X,Fll.4> A 510 

ENO A 511- 


BLOCK DATA B 1 

COMMON /JANl/ AJ,PIE,G,ZZ,PSTD,TSTD,PSIR,RA,DELO,SQTHl B 2 

DATA AJ/777. 97/, P I E/ 3 - 14159/ , G/ 32 . 17/, III. 0107663/, PSTD/2 1 1 6. / , T ST B 3 

1D/519./.PSIR/1. 5785/, RA/53. 331/ B 4 

END B 5- 


5 

6 
7 


11 

12 

13 


SUBROUTINE TURBIN ( T4 , P4PS0 , 0H54, AKC, AKU ,WC , ALPH AD, EN, F4, F I RST, D1 , D 1 
1U,DMT01,T5,P5P4,P5PS0,VMX5,HTR5, AMB AR , D5D1 , ETABAR, VX , VUl , VU2 , BET A1 0 2 


2 ,HTR4, 04D1 , VMX4 ,WM4, VM4, VXU ) D 3 

COMMON /JANl/ AJ,'PIE,G,ZZ,PSTO»TSTD,PSIR»RA,DELOfSQTHl D 4 

AMBAR=U*U/G/AJ/0H54 D 5 

ALPHA=ALPHA0/57 . 3 D 6 

C 0 7 

C CALCULATE VELOCITY DIAGRAM D 8 
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SUBROUTINE COMPRS ( T 1 , P2P1 , PIPSO, ETACO, HI , SWC1 , T2, P2PS0, DH2 1 , H2 ) C 1 

COMMON /JANl/ AJ,PIE,G,ZZ,PSTO,TSTO, PSI R ,RA, OELO , SQTH1 C 2 

C IF ETACO GT. 0 ETACO IS TREATED AS AN ADIABATIC EFFICIENCY C 3 

C IF ETACO LT.O ETACO IS TREATED AS A POLYTROPIC EFFICIENCY C 4 

P2PS0=P1PS0*P2PI C 5 

PH11=P0LY(T1,2) C 6 

IF (ETACO) 4,4,2 C 7 

2 PHI2I=PHI 1+AL0GIP2P1 )*RA/AJ C 8 

T2I=P0LY( PHI2I.8) C 9 

H2I=P0LY( T2I , 1 ) C 10 

DH21=(H2I-H1 l/ETACO C 11 

H2=H1+DH21 C 12 

T2=P0L Y(H2,7) C 13 

GO TO 6 c 14 

4 PHI 2 = PH 1 1 — ALOG ( P2P1 ) /ETACO*R A/ AJ C 15 

T2=P0LY(PHI2,8> ’ C 16 

H2=P0LY(T2,l ) C 17 

DH2l=H2-Hl C 18 

6 RETURN C 19 

END C 20- 


28 



c A 70 

IF < JCGEOM.EQ. 3) GO TC 14 A 71 

10 DHDTC2=I. /SORT ( 1 .+W1 AC1/ (TEM9*0HDTC1**2 ) I A 72 

DC2DCI=DHDTC1/DHDTC2 A 73 

GO TO 16 A 74 

12 TEM10=( (1.40HOTCD/2. 1**2 A 75 

TQ=Hl»;i/ (4.*TEH9*TEM10) A 76 

OHOTC2 =(l.-TQ)/(l.+TQ) A 77 

0C2DCI = ( 1 .+DHDTC1 ) /( 1 .+DHDTC2 ) A 78 

. GO TO 16 , ... A 79 

14 0HDTC2=SQRT< 1.-W1AC1/TEM9) ' , - ' A 80 

0C2DCI =1 . A 81 

16 DELB2=DELC2*PB2PB1 A 82 

DELTl=DELB2 A 83 

HT1 =POL Y ( TT1 , 1 ) A 84 

PSIHT1=P0LY(TT1,4I A 85 

HO=HF*- 1254. A 86 

F8=<HU-HC2)/(H0-(HC2/ETAB>-PSIHT1| A 87 

FAR=F3/(ETAB*(1.-FB/ETABI) A 88 

C A 89 

C HTIG ENTHALPY OUT OF COMBUSTOR A 90 

C CHRG8 CHARGEABLE BLEED A 91 

C HTIGT EFFECTIVE ENTHALPY AT ROTOR INLET A 92 

C HN EFFECTIVE ENTHALPY AT TURBINE EXHAUST AFTER COOLING AIR IS MIXEO A 93 

C DELHT TURBINE ENTHALPY DROP DUE TO WORK REQUIRED BY COMPRESSOR A 94 

C AND PUMPING OF COOLING AIR A 95 

C BUSER = BLEED FLOW! LBS/SEC ) /W l A 96 

C IF KIND= 0 PROGRAMMER FURNISHES COOLING BLEED A 97 

C IF KIND = 50 LIFT ENGINE A 98 

C IF KIND = 100 CRUISE ENGINE A 99 

C A 100 

INDEX=KIND/50 A 101 

IF (INDEX-1) 18,20,22 A 102 

18 PCA=PCA A 103 

GO TO 24 A 104 

20 PCA=.000110*TTl-.242 A 105 

GO TO 24 A 106 

22 PCA-.0001504TTl-.297 A 107 

24 CHRG8=.50 A 108 

8ETC0L=PCA*( 1. -BUSER) A 109 

BETTOT=BETCOL+BUSER A 110 

TEM11=(1.+FAR)*(1.-BETT0T) a 111 

htig*hti+fb*psihti a 112 

DELHB=HTI G-HC2 A 113 

UTM=0TM0CI*UTIPCC*SQTH1 A 114 

dhpump=utm*utm/aj/g a 115 

0ELHT=(0ELHC+DHPUMP*CHRGB*BETC0L)/(TEM11+( 1,-CHRGB)*BETC0L ) A 116 

hcolp=hcool*ohpump a 117 

26 W1=WG/(TEM11*BETC0L) A 118 

C A 119 

WT=WG A 120 

WB1=W1*< l.-BETTOT) A 121 

WB2=WB1»( l.+FAR) A 122 

AC 1=W1 *SQTH1 / < 0ELC1*SWC 1 ) A 123 

0CI=SQRT(4.*AC1/PIE> A 124 

DTM=0C1*0TM0C1 A 125 

C A 126 

C ROTOR INLET CONDITIONS A 127 

C A 128 

HTIGT=(HTIG*TEM11 + (1.-CHRG8)*BETC0L*HC00L)/(TEM1H-1 l.-CHRGB ) *BETCO A 129 

ID A 130 

FARR*FAR/( 1 .+8ETC0L* ( l.-CHRGB)/! l.-BETTOT ) ) A 131 

FR=ETAB*FARR/(l.<-FARR) A 132 

TTR=TT 1 A 133 

28 DEL TAT = (HTIGT-POLY(TTR, 1 )-FR*POL Y( TTR, 4 ) ) / ( POL Y ( TTR,- 1 ) + FR*POL Y( T T A 134 

1R.-4)) A 135 

TTR=TTR+OELTAT A 136 

IF (ABS(OELTAT).LT..l) GO TO 30 A 137 

GO 10 28 A 138 

30 WR=W1*(TEM11+BETC0L* I l.-CHRGB)) A 139 
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Figure 3. - Continued. 
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TS4=T4— V4*V4/ (2.*G*A J*CP4) 

D 

80 


VM4=V4/SQRT(GAMMA4*G*R4*TS4) 

D 

81 

53 

IF (FIRST. 60.0. ) GO TO 14 

D 

82 


VMX4=VX/SQRT <GAMMA4*G*R4*TS4) 

D 

83 

57 

GF2 =0. 5* ( GAMMA4+ 1 . 0 ) / ( GAMMA4— 1.0) 

D 

84 


W0N=(1.0+(GAMMA4-1.0) *0 . 5*VM4*VM4 ) **GF2 

D 

85 

58 

TW0=S3RT( GAMMA4*G/R4 ) 

D 

86 

59 

P4=P4PS0*PSTD 

0 

87 


AF4A1=4.*WC*SCRT( T4) *WON/ ( P4 *TWO*P I E*D1 *D1 *VMX4 > 

D 

88 

60 

RK=4. 0*DMTD1*DMTD1 

D 

89 


RK=AF4 Al/RK 

D 

90 


HTR4=( l. 0— RK ) / ( l.O+RK) 

D 

91 


D40 1=2. 0*DM TD1 / ( 1.0+HTR4) 

D 

92 


CP5=P3LY(T5,3)+F4*P0LY(T5,6) 

D 

93 

62 

P5PS0=P4PS0*P5P4 

D 

94 


GAMMA5=1.0/( 1.0-R4/AJ/CP5) 

D 

95 


TS5=T5-V5*V5/(2.*G*AJ*CP5) 

D 

96 


VMX5=VX/SQRT(GAMMA5*G*R4*TS5) 

D 

97 

64 

WM4=V5 /SORT ( G AMMA5*G*R4* TS4 ) 

D 

98 

65 

VM5=V5/SQRT(GAMMA5*G*R4*TS5) 

D 

99 

66 

GF2=0.5*(GAMMA5*1.0) / (GAMMA5-1.0) 

D 

100 


WON= ( 1 . 0* ( GAMMA5— 1 .0 ) *0. 5*VM5*VM5 ) **GF2 

D 

101 

67 

TWO=SORT ( GAMMA5*G/R4 ) 

D 

102 

68 

P5=P5PS0*PSTD 

D 

103 


AF5Al=4.*WC* SORT ( T5 ) *WON/ ( P5 *TWO*P I E*Dl*Dl*VMX 5 ) 

D 

104 

69 

RK=4.0*DMTD1*DMTD1 

D 

105 


RK= AF5 A 1 /RK 

D 

106 


HTR5 = ( 1 . O-RK ) / ( 1 . O+RK ) 

D 

107 


050 l=2.0*DMTDl/( 1.0*HTR5 ) 

D 

108 


RETURN 

D 

109 


END 

D 

110- 



FUNCTION POLY (X,N) E 1 

DIMENSION ICC(8), I0R0(8), CI300) E 2 

DATA ( ICC( I ) , 1=1, 8)/133, 146, 161,174,187,200,209,226/ E 3 

DATA ( IORDI I ) , 1 = 1, 81/4, 5, 4, A, 4, 2, 6, 6/ E A 

DATA (C(I ), 1=133, 2421/100. ,107. 2616, 2. 4975E-1, -2. 2658E-5 , l . 96075E- E 5 
18.-3.675E-12.2000. , 116.484, . 2096 10 , 2. 554713E-5 , -3. 338588E-9, 1.8843 E 6 
26-13,6 000. ,100. ,1.14 506, 1. 49413E-3 , - 1 .7983 IE-6, 1 .39476E-9 ,-5. 851 4E E 7 
3-13, 1.0156E-16, 1400. , 1 . 429 38 , 4.437 1 22E-4,- 1 . 4891 8E-7, 3.4627E- 1 1, -4 E 8 
4.513E- 15,2.476-19, 6000. , 100., . 25232, — 5. 44152 E- 5, 7. 0682E— 8 ,— 2 . 017 IE E 9 
5-ll,-5.1E-16,1400. ,. 1861 ,8. 01486-5,-2. 32 78E-8, 3. 41 635E-1 2, - 1. 989E- E 10 
616,6000., 100. ,9 57. 02 8, - 1.3972476-2,2. 728 088E-4, -7. 874997E-8, 1.13 11 E 11 
784E-1 1 ,2000. ,989.0299,2.4210586-2, 1.7198626-4,-2. 1315 IE-8, 1.027056 E 12 
86-12,6000. , l00.,-.224336,4.84768E-4,— 7.1 84546E-8, -4.31 5008E-12, 3.1 6 13 

9 347 44E -15, 2000. , -4. 089761,5. 5046576-3,-2. 446699E— 6, 4. 89935 16-10,-3 E 14 
$.6173836—14, 6000. , 100 . , 3 .986 078 E— 2 , 3 .562965E-4 , -6 • 5 1 2821 E- 8, 2000. , E 15 
$.2063979, 1.7643436-4,-1. 6831 37E-8, 6000., 100., -62. 6855 16, -6. 9221 196 E 16 
$E— 1,1. 9220060E— 2,-1. 1 34561 IE-5,-9. 8360870E-8 , 2 . 2458 196E- 10, - 1 .4509 E 17 
$(f78E-l3,614., 1946 .1462,— 8. 5195086, 2. 6225174E-2, -2. 6 54375 9E— 5,1.147 E 18 
$9667E-8,-6.42755396-13,-6.082041E-16,2000. ,1.4,415.37102,58.865753 E 19 
$,-312. 50006,16.528755,-29.769562,5.6792908*53.206712, 1.93384,553.8 E 20 
$4976,123.43567,-174. 32 12 9, 22 . 1 84326 ,-104. 70823, -43. 322374, 82.93804 E 21 
$9,2.3/ E 22 


DATA SOURCE — HALL AND WEBER, NACA RM E56B27 





POLY 

FUNCTIONS 



FOR 

AIR 






POLYI TEMP, 

1 > 

= 

ENTHALPY 

POLY (TEMP, 

4) 

= PSI (ENTHALPY! 

POLY(TEMP, 

21 

= 

ENTROPY 

POLYI TEMP, 

51 

= PSI (ENTROPY) 

POLY ( TEMP , 

3) 

= 

CP 

POLYI TEMP, 

6) 

= PSI ( CP ) 

POLYI TEMP, 

-1 ) 

= 

CP 

POLY ( TEMP, - 

4) 

= PSI (CP ) 

POL Y ( TEMP , 

-2) 

= 

CP/TEMP 




POLYIENTHALPY, 

7) 

= 

TEMP 

WHERE PSI IS 

AN 

INTERPOLATION 

POLY ( 5 NTROPY , 

8) 

= 

TEMP 

FACTOR (SEE 

RM 

E56B27 ) 


Figure 3. - Continued. 
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c 


J=I ABS ( M ) 

K= I CC { J ) 

L= I 0R3 ( J ) +2 

2 IF IX-CIK)) 4,6,6 

4 K=K-L 

GO TO 2 
6 K=K+L 

IF (X-C(KI) 8,8,6 
8 IF (M) 14,10,10 

10 L=L-l 

POL Y=C ( K— 1 ) 

DO 12 N=2,L 
KN=K-N 

12 POLY=POLY*X+C(KN) 

GO TO 18 
14 L=L-2 

POLY=FLOAT(L)*C«K-l) 

00 16 N=2,L 
KN= K-N 
LN=L-N+1 

16 POLY=POLY*X+FLOAT(LN)»C(KNI 
18 RETURN 

END 


E 23 
E 24 
E 25 
E 26 
E 2 7 
E 28 
E 29 
E 30 
E 31 
E 32 
E 33 
E 34 
E 35 
6 36 

E 37 
E 38 
E 39 
E 40 
E 41 
E 42 
E 43 
E 44 
E 45- 


2 


4 


6 


8 


10 

12 

C 

14 

16 


FUNCTION VISC IT) F 

DIMENSION A I 50 ) , V(50) F 

DATA ( VI J) ,J = 1, 50 1/7 3.8, 136. 0,185. 2, 227. 2, 264. 7, 299. 2, 33 1.3, 361.4, F 

1389.8.417.1.443.5.469.5.495.1.519.7.543.6.567.0. 589.8.612.1.633.9, F 

2655.3. 676.3.697.0. 717.3.737.3.757.0.776.5.795.6.814.5.833.2.851.6, F 
3869.8,887.8,905.6,923.2,940.6,957.9,974.9,991.8,1008.6,1025.2,1041 F 

4.6.1058.0. 1074.1.1090.2.1106.1.1121.9.1137.5. 1153. 1.1168.5. 1183.8, F 


TKELVN=T*5./9. F 
A ( 1 ) =1 00. F 
DO 2 N =2 , 50 F 
A ( N ) =A ( N-l 1 + 100. F 
CONTINUE F 
IF ( TKELVN. LE.A< 1 ) ) GO TO 8 F 
IF ITKELVN.GE.AI50)) GO TO 10 F 
K=2 F 
IF (TXELVN.LT. A ( K ) ) GO TO 6 F 
K=K+ 1 F 
GO TO 4 F 
TOP=A (K)-TKELVN F 
DELTT=T0P/100. F 
DV= V ( X )— V (K-l ) F 
01 FF=OELTT#DV F 
V I SM=V ( K ) — DI FF F 
V I SC=V ISM/ 14882000. F 
GO TO 12 F 
VISC=V(1)/14882000. F 
WRITE (6,14) F 
GO TO 12 F 
VISC=V (50)/14882000. F 
WRITE (6,16) F 
RETURN F 

F 

FORMAT (32H TKELVN IS LESS THAN 100 DEGREES) F 
FORMAT ( 36H TKELVN IS GREATER THAN 5000 DEGREES) F 
END F 
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Figure 3. - Continued. 
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SUBROUTINE DUCT2 ( PI , P2P 1 , AMD2 , TO, WOW 1 , EQWC l , DELC1 , GAMMD , RD , P2ST AT 6 1 


I.AD2A1) G 2 

COMMON / JAN1/ AJ,PIE,G,ZZ,PST0,TSTD,PSIR,RA,DEL0,SQTH1 G 3 

P2=P1*P2P1 G 4 

TEMl = l.-MGAMMD-l. )/2.*AMD2**2 G 5 

P2STAT=P2/(TEMl)**(GAMMD/( G AMMO- 1.1) G 6 

SRTDTF=SQRT(TO)/(SQTHl*22. 78166) G 7 

TEM2=W0W1*EQWC1*SRTDTF/P2*DELC1*PSTD G 8 

TEM3=AMD2*SQRT(GAMMD*G/RD) G 9 

TEM4=TEM1**( (GAMMD+l. ) / < 2. *( GAMMD- 1 . > ) 1/TEM3 G 10 

AD2 A 1 = TEM2*TEM4*ZZ G 11 

RETURN G 12 

end G 13- 


Figure 3. - Concluded. 
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Figure 4. - Sample computer printout from gas generator/lift engine program. 
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E-6968 


DUCT/NOZZLE 0.9900 0. 02586 23.95 1229 




be aeronautical and space activities of the United States shdl be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shdl provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof." 

— National Aeronautics and Space Act of 1958 
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TECHNICAL TRANSLATIONS: 
published in a foreign language considered 
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SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 


TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
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contract or grant and considered an important 
contribution ro existing knowledge. 
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